arXiv:1501.01533v2 [math.AP] 19 Jan 2015 


The enclosure method for the anisotropic 
Maxwell system 

Rulin Kuan; Yi-Hsuan Lin; Mourad Sini 


Abstract 

We develop an enclosure-type reconstruction scheme to identify pene¬ 
trable and impenetrable obstacles in electromagnetic field with anisotropic 
medium in R 3 . The main difficulty in treating this problem lies in the fact 
that there are so far no complex geometrical optics solutions available for 
the Maxwell’s equation with anisotropic medium in R 3 . Instead, we derive 
and use another type of special solutions called oscillating-decaying solu¬ 
tions. To justify this scheme, we use Meyers’ L v estimate, for the Maxwell 
system, to compare the integrals coming from oscillating-decaying solu¬ 
tions and those from the reflected solutions. 

Keywords: enclosure method, reconstruction, oscillating-decaying solutions, Runge 
approximation property, Meyers L p estimates. 

1 Introduction and statement of the results 

Let SI be a bounded C °°-smooth domain in R 3 with connected complement 
R 3 \ S2 and flbea subset of O with Lipschitz boundary. We are concerned with 
the electromagnetic wave propagation in an anisotropic medium in R 3 with the 
electric permittivity e = (e^ (a;)) a 3 x 3 positive definite matrix and e(x) = eo( x) 
in f 1\D. We also assume that e{x) = eo(x) — £d(x)xd(x) with eo € C°°(Ci) a 
positive definite 3x3 symmetric matrix and e d (x) is a positive 3x3 symmetric 
matrix and /i a smooth scalar function defined on fl such that there exist /r c > 0 
and e c > 0 verifying 
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fi(x) > n c > 0 and E > e c |£| 2 W G R 3 , Vz G fl. (1.1) 

*-i=i 

If we denote by E and H the electric and the magnetic fields respectively, then 
the electromagnetic wave propagation by a penetrable obstacle problem reads 
as 

{ V x E — ik^iH = 0 in fi, 

V x H + ikeE = 0 in SI, (1.2) 

v x E = / on dn, 

with e = eo — £dXd, and the one by the impenetrable obstacle as 

{ V x E — ik^H = 0 
V x H + ikeE = 0 
vxE = f 
v x H = 0 


in fl\D, 
in f l\D, 
on <9f1, 
on dD, 


(1.3) 
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where v is the unit outer normal vector on dfl U dD and k > 0 is the wave 
number. In this paper, we assume that k is not an eigenvalue for CE3 and (113- 
Impedance Map: We define the impedance map A d : TH 2 (<9f 1) —>• TH 2 (<9fl) 
by 

A d (v x H |an) = {v x E\gn), 

where TH~i (dCl) := {/ € H~i(dft)\v ■ f = 0} and x is the standard cross 
product in R 3 . We denote by A@ the impedance map for the domain without 
an obstacle. 

Consider the anisotropic Maxwell system 

J V x E — ik[iH = 0 in £1, 

[V xH + ikeE = 0 in Q, ^ ' ’ 

where fi and e satisfy ED- We are interested in the question reconstruct¬ 
ing the shape of D using the impedance map Ap. This geometrical inverse 
problem is quite well studied in the literature see 0 and several methods have 
been proposed to solve it. In this paper, we focus on one of these method, 
called the enclosure method, which is initiated by Ikehata, see for examples 
[U [3] , and developed by many researchers [71 |U HU HU HU US] , [SI HU f° r the 
acoustic model, Eli for the Lame model and mm for the Maxwell model. 
The testing functions used in [UEU are complex geometric optics (CGO) solu¬ 
tions of the isotropic Maxwell’s equation. The construction of CGO solutions 
for isotropic inhomogeneous Maxwell’s equations is first proposed in [T7]. Af¬ 
ter that, the authors in [5] also constructed CGO solutions for some special 
anisotropic Maxwell’s equations. However, there are not yet of CGO solutions 
for general anisotropic Maxwell system. Besides, CGO solutions, another kind 
of special solutions for anisotropic elliptic system was proposed for substitution 
in m and m • They are called oscillating-decaying (OD) solutions. Inspired 
by HU and EE our idea is to reduce ed to an elliptic systems and then 
use the results in |T5] to construct oscillating-decaying type solutions to the 
anisotropic Maxwell system. Precisely, we can decompose the equation ED 
into two decoupled strongly elliptic systems. The main difference between the 
construction of the oscillating-decaying solutions in m and ours is about the 
higher derivatives of oscillating-decaying solutions. 

One of the main differences between the CGOs and the oscillating-decaying 
solutions is that, roughly speaking, given a hyperplane, an oscillating decaying 
solution is oscillating very rapidly along this plane and decaying exponentially in 
the direction transversely to the same plane. Oscillating-decaying solutions are 
special solutions with the phase function having nonnegative imaginary part. In 
addition, these oscillating decaying solutions are only defined on a half plane. 
To use them as inputs for our detection algorithm, we need to extend them 
to the whole domain O. One way to do the extension is to use the Runge 
approximation property for the anisotropic Maxwell’s equation. The Runge ap¬ 
proximation property will help us to find a sequence of approximated solutions 
which are defined on fl, satisfy ED and their limit is the oscillating-decaying 
solution. Note that it was first recognized by Lax m that the Runge approxi¬ 
mation property is a consequence of the weak unique continuation property. In 
EE the authors already proved the unique continuation property and based on 
it we derive the Runge approximation property for the anisotropic Maxwell’s 
equation. 
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To be more precise, let w be a unit vector in R 3 , denote f2 t (w) = Q,C\{x\x-vj > 
t}, E t (w) = n fl {x\x ■ u> = t} and set ( E t ,H t ) to be the oscillating-decaying 
solution for the anisotropic Maxwell’s equation in 

Support function: For p S § 2 , we define the support function of D by hr>(p) = 
inf xeD x ■ p. 

When t = ho(p), which means E t (w) touches dD , we cannot apply the 
Runge approximation property to ( E tj H t ) in D t (w). Therefore, we need to 
enlarge the domain D t (w) such that the OD solutions exist and the Runge 
approximation property works. Let p be a positive real number, denote 
and T, t -n(oj) and note that Li t _ v (uj) C D t (w) Mp > 0. We can find (E t _ v , H t _ v ) 
to be the OD solution in J7 t _ 7J (o;). By the Runge approximation property, there 
exists a sequence of functions t, H Vt t)} satisfying the Maxwell system in Q 
such that (E V: e,H Vt e) converges to (E t _ v ,H t _n) as t —> 00 in L 2 (f7 t _ I) (w)) and 
in H(curl,D) by interior estimates since D <s Q t - V (uj). In addition we show 
that (Et-rjjHt-rj) converges to ( E t ,H t ) in H(curl,D) as p —> 0. Then we can 
define the indicator function as follows. 

Indicator function: For p G S 2 , t > 0 and t > 0 we define the indicator 
function 

/ p (r,f) := lim lim I n /(T,t), 

r )—¥0 l^too H 

where 

I'I’ £ {r,t) := ikr [ {v x H v ^) ■ {(h D - A 0 )(V x H v> i) x v)dS. 

Jd n 

Goal: We want to characterize the convex hull of the obstacle D from the 
impedance map A£>. 

The answer to this goal is the following theorem. 

Theorem 1.1. Let p G S 2 . For the penetrable (or impenetrable) obstacle case, 
we have the following characterization ofhu(p). 

j linv^oo \I p (r,t)\ = 0 when t < ho{p), 

|liminf T _ ) . 00 \I p (r, h D (p))\ > 0, 

To prove Theorem I f .11 for the penetrable obstacle case, we need an appro¬ 
priate L p estimate of the corresponding reflected solution. We follow the idea in 
[7] to prove a global L p estimate for the curl of the solutions of the anisotropic 
Maxwell’s equation, for p near 2 and p < 2. 

To prove Theorem ll.il in the impenetrable obstacle case, we use layer poten¬ 
tial arguments as in [7] coupled with appropriate L p estimates. Precisely, first, 
we use the well-posedness for an exterior isotropic Maxwell’s system with the 
Silver-Miillcr radiation condition and, in particular, the layer potential theory 
to find a suitable estimate for the solution of this exterior problem. Second, we 
decompose the reflected solution into two functions, one satisfies the reflected 
Maxwell’s equation with a zero boundary data, the other satisfies the original 
anisotropic Maxwell’s equation with the same boundary conditions which come 
from the reflected equation. For the first decomposed function, we use the L p 
estimates, and for the second function, we will use the wcll-posedness, in L 2 , 
for the anisotropic Maxwell’s system. Combining these two steps, we derive the 
full estimate for the reflected solution in the impenetrable obstacle case. 
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This paper is organized as follows. In the section 2, we give decompose the 
anisotropic Maxwell system into two strongly elliptic systems. In section 3, we 
use the elliptic systems derived in the section 2 to build the oscillating-decaying 
solutions for the Maxwell system. Then, we give the Runge approximation 
for the anisotropic Maxwell equation in section 4. In section 5, we prove the 
Theorem D for both penetrable and impenetrable obstacle case. Finally, in 
the last section, as an appendix, we provide some technical details which we 
postponed in the main text and recall some useful estimates for solutions of the 
Maxwell system. Before closing this introduction, let us mention that in the 
whole text whenever we use the word smooth it means C^-smooth. 


2 Reduction to strongly elliptic systems 

Our goal is to construct the oscillating-decaying (OD) solution for the following 
anisotropic time-harmonic Maxwell’s system 

V x E = ikpH 

V x H = —ikeE , , 

i div( e E)=0 ’ (21) 

div(nH) = 0, 

where E 1 H denote the electric and magnetic field intensity respectively, and g, 
denotes the positive scalar permeability, e denotes the permittivity, which is a 
real, symmetric, positive definite 3x3 matrix. 

Inspired by na, the first step of constructing OD solutions is to reduce m 
to a strongly elliptic system. In fact, we reduce the anisotropic Maxwell’s system 
m to two separate strongly elliptic equations (12.31) , while in na the isotropic 
Maxwell’s system is reduced to an elliptic (a single Schrodinger) system with 
coupled zero-th order term. The following theorem is our reduction result. 

Theorem 2.1. We set E and H of the following forms 

E = ——e -1 V x (m _1 (V x B)) - e -1 (V x A) 

i k 

H = -/r^V x (e -1 (V x A)) - /^(V x B) 
k 

with A, B satisfying the strongly elliptic systems 

pVtr(M A \/A) - V x (e -1 (V x A)) + k 2 pA = 0 
eVtr(M B VR) - V x (^ _1 (V x B)) + k 2 eB = 0 

where M A ,M B are introduced in Theorem \2.4\ then E and H satisfy 

Remark 2.2. Theorem 2.1 shows that, if we can find solutions of (12.31) . then we 
can find solutions of m- 

Proof. In this proof, we will show the process of the reduction. And the proof 
that the systems (12.31) are strongly elliptic systems will be postponed to Theorem 

CO 


( 2 . 2 ) 


(2.3) 


4 




As in ns* we set the following two auxiliary functions which are similar to 
what they used: 

$ = -j-div(e-E’) 
k 

and 

'L = ydiv(nH). 
k 

Note that $ and T are actually zero by the Maxwell’s equation. We consider 
the following first-order matrix differential operator P 


P = 


Note that P is a 8 x 


( °1 
Ji _1 v 

div(e(-)) 

0 

0 

Vx 

° ^ 
0 

0 

-Vx 


0 

e -1 V 

V o 

0 

div(M-)) 

o / 

8 matrix. 

Let 





Y = 

(*\ 

E 

H 





V W 




Then the problem (12.11) can be rewritten as follows: 


PY = —ikVY, 


where 

/ 1 0 0 0 \ 

0 e 0 0 

V ~ 0 0 n 0 

\ 0 0 0 1 / 

Thus, the Maxwell’s system m implies 

(P + ikV)Y = 0 and $ = ^ = 0. (2.4) 


It is easy to see that conversely (G3D implies the Maxwell’s system, and hence 
they are equivalent. 

The first idea of the reducing process is to construct a suitable Q , which can 
make (P + ikV)Q a “good” second-order differential operator. Then, a solution 
X for the problem 

(P + ikV)QX = 0 (2.5) 


will give rise to a solution Y = QX for 

{P + ikV)Y = 0. 

Moreover, if we find the solution X such that the first and the last component 
of Y = QX are zero, then we obtain solutions for the Maxwell’s system. 
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We try the matrix differential operator Q = Q — ikl, where 

/ o 

V 

0 

V o 


Q = 


div(e(.)) 

0 

-0 _1 (V x (•)) 

0 


0 0 \ 

^(Vx(')) 0 

0 V 

div(/x(-)) 0 ) 


( 2 . 6 ) 


Then 


(.P + ikV)Q 
= (P + ikV)(Q-ikI) 

= PQ — ikP + ikVQ + k 2 V 
( div(eV) 0 0 0 


0 

0 

0 

0 


L\ 

0 

0 


0 

L2 

0 


—ikyL 

0 

0 


V 

/ 0 z/cdiv(e(-)) 


div(/iV) / 

—ifcdiv(e(-)) 0 

0 -ifcV x 

ikV x 0 

0 — ikdiv(n(-)) 

0 0 \ 


0 

0 

—ifce _1 V 

0 


ikeV 

0 



ikS/ x 

0 

0 

— ikS7 x 



0 

ikfiS7 

0 

0 


ikdiv(n(-)) 

0 

k 2 

0 

0 

0 




0 

k 2 e 

0 

0 



0 

0 

fc 2 /i 

0 




0 

0 

0 

k 2 

J 




V 

( div(eV) + k 2 
ik(e — /r -1 )V 
0 
0 


0 

Li + k 2 e 
0 
0 


0 

0 

L 2 + k 2 ^i 
0 


\ 


ik(n — e 1 )V 
div(/iiV) + k 2 J 


where 


Li = At - 1 V(div(e(-))) - V x (/i -1 (V x (•))) 
L 2 = e - 1 V(div(M(-))) - V x (e -1 (V x (•))). 


(2.7) 

( 2 . 8 ) 


A prominent feature of the above operator is that it decomposes the original 
eight-component system into two four-component systems. Precisely, Set 


X = 


e 
h 

W 


then (12.51) can be separated into two systems: 

J div(eV<p) + k 2 ip = 0 

1 Lie + k 2 ee + ik(e — /r _1 )Vy> = 0. 
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and 


div(/iVi/0 + k 2 tp = 0 
L 2 h + k 2 f-ih + ik(/j, — e -1 )Vf/’ = 0. 

Moreover, 

Y = QX 

' / 0 div(e(-)) 

V 0 

0 V(VxB) 

\ 0 0 

/ div(ee) — ihip 

Vip + e _1 (V x h) — ifce 
—^ _1 (V x e) + V0 — ikh 
y div(/z/i) — iki/i ) 

Therefore, the problem of finding the solutions X of 

(P + ikV)QX = 0 with the first and last component of QX being 0 (2.9) 

is equivalent to the problem of finding solutions of the following two separate 
systems: 

div(ee) — iktp = 0, 

div(eV<^) + k 2 tp = 0, (2-10) 

p, _1 V(div(ee)) — V x (/i -1 (V x e)) + k 2 ee + ik(e — = 0, 


0 

e -1 (V x (•)) 
0 

di v(/x(-)) 

\ 


0 

0 

V 

0 


— ikl 


X 


and 


di v(iih) — ikip = 0, 

div(/uVi/>) + fc 2 i/i = 0, (2-11) 

e^ 1 V(di y{fih)) — V x (e^ 1 (V x h)) + k 2 fih + ik(fj, — e^ 1 ) Vi/; = 0. 

Notice that if we set e in the following form 


e=--(V ¥ ; + e- 1 (Vx^)), 


( 2 . 12 ) 


then the first equation of (12.101) becomes the same as the second one. For the 
third equation, we have 

/i _1 V(div(ee)) — V x (/r _1 (V x e)) + k 2 ee + ik(e — 


= -^ _1 V (div(eVip)) + jV x (/i" 1 [V x (e _1 (V x A))] 

rC fC \ 

— ike\7<p — *fc(V x A) + ikeVip — -Y _1 V ( 'k 2 (pJ 


= V (div(eVyj) + k 2 ip) + -Vx (^t -1 [V x (e _1 (V x A))] ) - ikV x A 

= 0+jVx U _1 [V x (e _1 (V x A))] ) - ikV x A, 
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by the second equation of (12.101) . Thus, by letting e be of the form (12.121) . the 
system (12.101) reduces to 


div(7Vy>) + k 2 (p = 0, 

V x ^/t- 1 [V x (e -1 (V x A))] - k 2 A^j = 0. 


(2.13) 


Similarly, by letting 


h = —-(V^i + /x -1 (V x B)) 


for some vector field B, we can reduce (|2.11l) to the following system: 
div(/rVV’) + k 2 if = 0, 

V x ^e -1 [V x (/x-^V x B))] - k 2 B^j = 0. 


(2.14) 


To resume, if we can find solutions p,A,if and B of (12.131) and (12.141) . we can 
find solutions of the problem (12.91) and therefore the original problem (12.11) . 

Now let us focus on (12.131) and (12.141) . The goal is to find special solutions 
(e.g. oscillating-decaying solutions) of (12.131) and (12.141) . The idea of doing that 
is to subtract zero terms of the form V x (Vlr(M^Vd)) and V x (Vfr(M s VI?)) 
from the second equations of (12.131) and (12.141) for some matrices M A ,M B , so 
that they become V x (C A A) = 0 and V x ( C B B ) = 0 with C A and C B being 
strongly elliptic operators. Precisely, we want to find suitable matrices M A and 
M b such that 


pVtr{M A VA) - V x (e _1 (V x A)) + k 2 pA = 0 (2.15) 

and 

eS7tr(M B S7B) - V x (^ _1 (V x B)) + k 2 eB = 0 (2.16) 

are strongly elliptic systems. In fact, by letting M A = m/i _1 i and M B = 
mpT 1 e, we can show that (12.151) and (12.161) are strong elliptic systems for arbi¬ 
trary positive constant to. The proof are given in Theorem 12.41 □ 

To prove Theorem l2.41 we start with the following computational lemma. 

Lemma 2.3. Let M be a matrix-valued, function with smooth entries and F be 
a vector field. Then the i-th component of the vector V x (M(V x F)) is given 
by 


(V x (Af(V x F))) i = Y, C i]k id je f k + R u (2.17) 

j,k,£ 


where 


Cijk@ — djphlki T dikMgj SjkMgi S,j Atj-j T ( 8, p, 6 ,jf ) tv (AT). 

and Ri contains the lower order terms. Here, Sij is the Kronecker delta, Mij is 
the ij-th entry of M, and F = (/i, / 2 , /a) T . 






























Proof. We prove it by direct computations. For any vectors a, b, letting c = 
a x b, we have 


Cm — ^ ^ 
ki 

where a = (ai,a 2 ,a 3 ) T , b = ( 61 , 62 , & 3 ) T , c = (ci,C 2 ,C 3 ) t and £ m tk denotes the 
Levi-Civita symbol. Therefore, we obtain the m-th component of V x F: 

(VxFj =^mMk. 

' ' m kt 

Then, the n-th component of M(V x F) is 


(m(V x F) j = M nrn £ m e k defk■ 

\ / n ™ u p 


n m,k,l 

Finally, taking the curl operator on the vector M(V x F), the i -th component 
of the resulted vector is 

V x (M(V x F)) .= E £ijn&j (MrimSrnlkdpfk) 

/ 1 j,n,m,k,£ 

= ^ ^ £ijn£m£k((&jMnrnjdpfk T 


j,n,m,k,£ 


Thus 


V x (M(V x F)) ). = E Cijkidjifk T Ri 5 

' i j,k,£ 


where 


Cijkl •— ^ ^ £jjn£m£kM n mi Ri •— ^ ^ ^ijn^m£k{^jM n m)dzfk- 

m,n j,m,n,k,£ 


Since 


^ijn^-mik — 


&im $ik 

3jm $jl 3jk 

&nm ^nt $nk 


— ^im\^j£^nk $ni$jk) &i£ jm$nk $nm$jk) T $ik \^jm^n£ 5 

we can obtain 

Cijk£ = ^ ^ f $im{$j£$nk $n£$jk) $il{$jm$nk $nm$jk) 


T $ik(fijm$n£ &nm &jZ) J M n m 

= ( SjeM ki - SjkMa) - SuMkj + 5tf6jktr(M) + 5 ik Mej - 5 ik 5jetr(M) 

= A4 kl T 6zi;A/ f-j 5j k Mpi djjAI k j T (S l( s :l k &ik ) tr ( A/) . 


□ 
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Theorem 2.4. Assume that p, is a smooth, positive scalar function and e is 
a symmetric, positive definite matrix-valued function with smooth entries. The 
eigenvalues of e are denoted by Ai(x), A 2 (x) and \z(x). Assume there exist 
positive constants p.o, A, A such that for all x £ 

0 < n(x) < fj, 0 


0 < A < Ai(x) < A 2 (x) < Aa(x) < A. (2-18) 

Then S2.15 \) and \2.16\) are uniformly strongly elliptic by letting M A = m/W 1 / 
and M B = for arbitrary positive constant m. Here I denotes the 3x3 

identity matrix. 

Proof. To see whether (12.151) and (12.161) are strongly elliptic, we only have to 
check the leading order terms of (12.151) and (12.161) . We divide this proof into 
two parts, Part A and Part B, to deal with the equation (12.151) for A and the 
equation (12.161) for B respectively. 

Part A. By Lemma [2731 

(^iS7tr(M A \7 A) - V x (7 _1 (V x A))^ 

= £ /aSijdj ( M A d e A k ) ~ Y. - R? 

jkt jkt 

= Y {^M A - C A kt )dpA k + Y d 3 M A )d,A k - R a 

jkt jkt 

= Y C ^ed 3 eA k + Y liSiMMfadeA, - R A , 

jkt jkt 

where C A ke = p,Sij M A k — C A kf are the coefficients of the leading order 
terms of (12.1511 and 

Cij k t = )ki T di k (e Sj k (e bi(.(e A (digdjk di k djz)tr(e ). 

Recall that (12.151) is called uniformly strongly elliptic in some domain £l 
if there exists a positive cq > 0 independent of x £ fl such that 

Y C Hkti.x)a l a k b j b ( . > c 0 |a| 2 |b| 2 (2.19) 

ijkt 
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for any a,bgl 3 and for all x € fl. Now 

F Cijke a i a kbjbe = ^2 (V'SijMzk ~ &’ijkt ) CLiCLkbj bp 
ij ki ij kl 

= //(a • b)(b T M A a) 

^ ^ )ki ~ I - bik(t )lj bjk (c 

ijki ^ 

)&? H“ {biZ^jk ^ik^ji)bv{c did^bjbz 
= ^(a-b)(b T M A a) 

- ^|b| 2 (a T e _1 a) + |a| 2 (b^e” 1 b) - (a ■ b)(b T e - 1 a) 

— (a ■ b)(a T e _1 b) + fr(e _1 )(a ■ b) 2 — fr(e _1 )|a| 2 |b| 2 J 

= fr(e^ 1 )|a| 2 |b| 2 — |a| 2 (b T e _1 b) - |b| 2 (a T e -1 a) - tr(e~ 1 )( a • b) 2 
+ 2 (a • b) (b T e^ 1 a) + /i(a • b) (b T M 4 a) 

since e (and hence e -1 ) is symmetric. Let S be the orthogonal matrix 
such that e = S T DS 1 where D = diag(Ai, A 2 , A 3 ). Thus e _1 = S T D~ 1 S. 

Also let M A = S T N A S. By letting v = S'a/lal and w = 5 b/|b|, it’s easy 
to see that (12.191) holds for all a, b £ R 3 iff 

fr(e -1 ) - (w r D~ 1 w) - (v T D~ 1 v) - tr (e" 1 )^ • w ) 2 
+ 2 (v • w) (w T F - 1 v) + /i(v • w) (w T W 4 v) > co 

for all v, w € R 3 such that |v| = |w| = 1. Note that fr(e _1 ) = tr(H _1 ) = 

A )” 1 + A ^ 1 + AJ 1 . In summary, we find that (12.151) is uniformly strongly 
elliptic on Q iff 

inf f min F(v, w)| >0, (2.20) 

xeo Vlv|=|w|=i y 

where 

F(v, w) = ^tr(D~ 1 ) - (w t F _ 1 w) - (v T F>- 1 v) - tr(D _1 )(v • w) 2 

+ 2(v ■ w)(w T Zl _1 v)^ + /x(v • w)(w t N a v') 

=: G(v, w) + //(v ■ w) (w T IV' 4 v) . 

We will show that 

G(v,w) > Aj 4 (l - ( V • w ) 2 ) (2.21) 

under the constraints |v| = |w| = 1. Then, by choosing M A = m^” 1 / for 
some positive constant m, we also have N A = and 

F(v, w) = G(v, w) + m(v ■ w) 2 

> A 3 1 (1 - (v • w) 2 ) + m(v • w) 2 
= A3 1 + (m — A3 4 )(v • w) 2 . 
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Now since 0 < (v • w) 2 < 1, if m > AJ 1 , we have F(v, w) > A^" 1 , while 
if m < A^" 1 , we have F(v, w) > A^ -1 + (m — A^" 1 ) = mn. Remember that 
A^ 1 (a;) > A -1 on f2, we conclude that F(v,w) > min(A _1 ,m) for all 
|v| = |w| = 1 and all x G fi. 


It remains to show (12.211) . For this, note that 

G(v, w)= ^ E i l ~ w j ~ v ) ~ ( v ' w ) 2 + 2 ( v ' w H™i) = : E X J lR j- 

j= 1,2,3 ' ' j 

We can prove Kj > 0 as follows: Since (v ■ w) — v\W\ = V 2 W 2 + V 3 W 3 , by 
Schwarz inequality we have 


|(v • w) — v\W\\ < \jvl + vl^jwl + wl = -wf. 

Taking square, we obtain 

(v • w) 2 — 2 (v • w)uiu;i + v\w\ < 1 — vf — wf + v\iu\, 

which means K\ > 0. Similarly K 2 ,K 3 >0. As a consequence, since 
Ar 1 > A(T 1 > AJ" , we have 

G(v, w) > A 3 1 (ATi + K 2 + K 3 ) = A ^" 1 (l - (v • w) 2 ), 

which completes the proof of Part A. 

Part B . For (12.161) . we have 


( 7 Vfr(M B VB) - V x (m _1 (V x B)) 


- E ( M ?kdeB k ) - E C? m d jt B k - Rf 

jki jki 

= E ~ Cf jkt )d jt B k + E 7 - Rf 

j ki j ki 


( 2 . 22 ) 


where 


^ijki — 8jpfl 8 k i A fiikfJ* $tj ( \jk 7 8gi 

8 it\A 8 k j A ( 8 t p 8 jk 8 i k 8 jfj ^k(jl /) 


= M 


-1 




Denote the coefficients of the leading order terms of (12.221) by Cf M , we 
have 


°ijki 


Mtk ~ Cijkt — e ijM ek — /j ( 8 uSjk 



By choosing M B = m/i x e we obtain 

E C ij k eaiakbjbe = iT 1 fm(a T 7b ) 2 - ((a • b ) 2 - |a| 2 |b| 2 ) 

ij ki ^ 
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for all a, b £ R 3 . Remember that e = S T DS. Since we have assumed 
/r -1 > fi o for some positive constant fj, o, by letting v = .Sa/fyl and w = 
Sb/|b| for a,b/0, we see to prove C B ki aiakbjbe > co|a| 2 |b| 2 for some 
constant cq > 0 is equivalent to prove 

inf min iJ(v, w) > 0 , ( 2 . 23 ) 

x6fi|v| = |w| = l 

where U(v,w) = m(v T Dw ) 2 + (l — (v • w) 2 ) . Although (12.231) looks 
simpler than (12.201) . we fail to find a simple method as before to get a 
clear lower bound. Nevertheless, it is also easy to see that (12.231) is true 
by continuity, as follows: If (v • w) 2 = 1, then v = ±w, and 

m(v T Dw ) 2 = m(Xivl + A2w| + A3U 2 ) 2 > mA 2 . 

By continuity, there exists s > 0 such that for 0 < 1 — (v • w) 2 < e 
we have m(v T Z?w) 2 > m\\/2. Thus for 0 < 1 — (v • w) 2 < e we have 
Lf(v, w) > m\\/2. While for 1 — (v • w) 2 > e, H(v, w) > e. Thus under 
the constraints |v| = |w| = 1 we obtain 

H(v, w) > min(?nA 2 /2, e) > min(mA 2 /2, e), 

where recall that A is the lower bound of Ai(x) on f l. This completes the 
proof of Part B. 


□ 

Remark 2.5. One can check that the C A and C B satisfy C A - k( . = G Uij and 
G Fjki = G kUj- And, by choosing M A = m/i -1 1 and M B = m/r _1 e as above, 
the C A and C B also satisfy such symmetry. This additional property is useful 
in the next section. 


3 Construction of oscillating-decaying solutions 

In this section, we will use the reduction results in section 2 to construct 
oscillating-decaying solutions of m- From now on, we suppose that /i > 0 
is a C°° scalar function and e is a 3 x 3 real positive definite matrix-valued 
smooth functions (i.e. every entry is a real C°° function) and E , H satisfy 

j V x E — ik^xH = 0 in fi, 

| V x H + ikeE = 0 in Q. 

In order to obtain the oscillating-decaying solutions of E and Eb, we have to 
construct the oscillating-decaying solutions for A and B. We follow the proof 
in m to construct the oscillating-decaying solutions for A and B , but here we 
need to derive higher derivatives for A and B. 

From he we borrow several notations as follows. Assume that C R 3 is 
an open set with smooth boundary and w £ S 2 is given. Let 77 £ S 2 and ( £ S 2 
be chosen so that { 77 , (, w} forms an orthonormal system of R 3 . We then denote 
x' = (x- 77 ,x-C). Let t £ R, flt(w) = I2n{x-u; > t} and E t (w) = fl D {x ■ w = f} 
be a non-empty open set. 
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Theorem 3.1. Given {rj,f,uj} an orthonormal system of R 3 , x' = (x ■ r], x • Cf) 
and f£l. We set fl t (w) = finjx-cu > t} and £ t (u;) = = t}, then We 

can construct two types OD solutions for the Maxwell system in Tl t (u>) which 
can be useful for penetrable and impenetrable obstacles respectively. There exist 
two solutions of iTOI) of the forms. The first one is 


' E = F\(x)e iTX <e- T ^- u -^ A t ( ~ x ">b + T 
H = Fl{x)e iTX ^e~ T( - x ' ul - t '> A ^ < - x ''>b + T 


Xt,b,t,N,u>( X i T ) + r X t,b,t,N,u>( X ’ T ) 


A,2 


Xt ,b,t,N,J X > T ) + r xub,t,N,J X ’ T ) 


in 

in f2t(w), 


where F\(x) = 0(t), F\(x) = 0(t 2 ) are some smooth functions and for \ a 
j = 1,2, we have 


(3.1) 
= 3, 


f \K:it,Nj X ^)h^ t{u) ) < CTMS/2 e -r(s-t) a ^ 

\ \\ r xub,t,N,w( x ’ r )llL 2 (n t (w)) < crJ- JV + 1 / 2 , 
for some positive constants oa and c. The second one has the form 

[E = + + r ’*£ ttN Jx, T ) 

\H = G 1 B (x)e iTX <e~ T ( x - u -V A ? ( - x ''>b + rf 


(3.2) 


Xt,b,t,N,J X ’ T ) + r xub,t,N,A X ’ T ) 


in f2t(a;), 
in n t (ui), 


where G B (x) = 0 {t),G 2 b {x) = 0(r 2 ) are some smooth functions and for \a 
j = 1 , 2 , we have 

j l|rf;f MiiV ,J^,r)|| i2(ntM) < cr l«|- 3 / 2 e -r 

\\\ r Xub,t,N,u( X ’ r )IU 2 (O t (w)) < CTl- N+1 / 2 , 


(3.3) 
= 3, 


(3.4) 


for some positive constants as and c. 

Proof. We want to find special solutions A, B € (C°°(12 t (w)\9£ t (w))nC 0 (f2 t (u;))) 3 
with r 1 satisfying Dirichlet boundary problems 


L a A -.= i_iS7tr(M A \7A) - V x (e'^V x A)) + k 2 p,A = 0 
A = e iTX < [xt{x')Q t {x')b + P A , t , b , N , u } 

and 

L b B := eVtr(M B VB) - V x (m _1 (V x B)) + k 2 eB = 0 
B = e ixx < { X t(x')Q t (x')b + P* tt>biNiU } 

where £ £ S 2 lying in the span of { 77 , <^} is chosen and fixed, Xt{x r ) £ C“(]R 2 ) 
with supp(xt) C £t(w), Qt{x') is a nonzero smooth function and 0 / £ C 3 and 

N is some large nature number. Moreover, f3 A t b t Nu (x', r), j3 B t b t N u {x', t) are 
smooth functions supported in supp(xt) satisfying: 

\\@Xt,b,t,N,uj{' 1 T ) ll-L 2 (R 2 ) < CT' 1 , ||/3f t , M ,jv,J-,T)|| L 2(R2) < CT -1 


in fli(aj) 
on £ t (w), 


(3.5) 


in fit(w) 
on £ t (u;), 


(3.6) 


for some constant c > 0. From now on, we use c to denote a general positive 
constant whose value may vary from line to line. As in |15l . A, B satisfy second 
order strongly elliptic equations, then it can be written as 


A — Ax t ,b,t,N,uJ 

B — 


= W X t ,b,t,N, u + r xt,b,t,N, u 
= W x t ,b,t,N, U + r xt,b,t,N,u: 
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with 


= Xt(x')Q t e^<e-^~'^') b + r^ t ^>, t ) 


(3.7) 


and r xtb,t,N^ r xtb,t,N,u satisfying 

\\ r x t ,b,t,N,J\HHn t (u)) < cr k ~ N+1/2 , \\rx tAt , N ,J H nn t (u> ) ) < cT k ~ N+1/2 , (3.8) 


where A A (-), Af(-) are smooth matrix functions with its real part ReA^(x') > 
0, ReAf (x') > 0 and h t jy , r x t ,b,t,N,u are a smooth functions supported 
in supp(xt) satisfying 


II d“ r xt ,b,t,At Jl £ 2 (SM^)) 

\\ d x r xt,b,t,N,J\^(n a ( u )) 


< c T l“l - 3 / 2 e -' r ( s ~ t ) a A 

< cr |a|-3/2 g -r(s— t)a B 


(3.9) 


for |a| £ N U {0} and s >t, where a A , clb > 0 are some constants depending on 
A a (x') and Af(x') respectively. We give details of the construction of A and 
B with the estimates (13.71) and (13.81) in the appendix. 

In Appendix 6.1, we derive the explicit representation of A and B. Recall 
that E and H are represented in terms of A and B as follows 


{ E = -t 7 _1 V x (/^(V x B)) - 7 -1 (V x A), 
H = ^ _1 V x (y-^V x A)) - m _ 1 (V x B). 


(3.10) 


Now, we can show that ( E , H ) satisfies (13.11) . (13.21) and we will use this form to 
prove Theorem 11.11 for the penetrable case. Similarly, we can show that (E, H ) 
satisfies (13.31) . (13.41) in order to prove Theorem ll.ll for the impenetrable case. All 
we need to do is to differentiate A and B term by term componentwisely. For 
the main terms of A and B, we can differentiate Xt(x')QtF ,TX '^e~ T ^ x ' u ~ t ^ At ( x )b 
and xt{x')Qte lTX ^e~ T( - x ' ul - t ^ A t ( - x '' > b directly and it is easy to see that 


V x A = TF A (x)e iTX ^e~ T( - x - ul - t ^ ( - x ^b + V x r A b t N Jx,T) + V 

V x B = rF B (x)e iTX <e- T( - x - ul - t ^ < - x ^b + V x r% t , b , t , N ,Jx,T) + V 


X r xt,b,t,N,v’ 
X r xt,b,t,N,ui > 


where F A {x) and Fb(x) are smooth matrix-valued functions and support in 
supp(xt(® / ))- F° r the penetrable obstacle case, we choose A = w A t btNu + 
r xt,b,t,N,ui be the oscillating-decaying solution satisfies L a A = 0 and B = 0 
(also satisfies L^O = 0) in fl t (u;), then (13.101) will become to 


E = — 7 _1 (V x A), 

i7=^- 1 Vx( 7 - 1 (VxA)), 


which means 


E = F\(x)e iTX <e- T ^ 
H = F\(x)e iTX '^e~ T ( : 


t) A ?{x')i) -p 

-t) A t( x ')b + 


r iHb,t,N, U ( X > T ) + r xub,t,N,J X ’ T )i 

T xXt,N,S x i T) + r^ 2 btNul {x, r), 
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where F\(x), F\(x ) are smooth functions consisting p(x), e(x), Q t (x '), Af(x') 
and their curls (it can be seen by directly calculation). Moreover, by suitable 
choice of b (for example, we can choose b ^ 0 is not parallel to £), we will get 
F\(x) = 0{t) and F\(x) = 0(r 2 ). Moreover, T^ b t Nu] and T^ b t Nu satisfy 

(ESI) for H = 1 and |a| = 2 , respectively, r^ btNuj and r^ btNu] satisfy 63 
for k = 1 and k = 2, respectively. Similarly, for the impenetrable obstacle case, 
we choose A = 0 and B = w^ t b t N ul + r^ t b t N ul in £2 t (w), then 

\e = G% (^yrx-f e -r(x-u,-t)Af (»') b + T^ N Jx, r) + rfj^jx, r), 

H = + + r^ N Jx,r), 


where G B (x) = 0{t) and G B (x) = 0(r 2 ) and 


■ Xtd 
r 2\ 


\a\=j and r^'\ 


Xt ,b,t,N,ui 


satisfies (13.91) for k = j. 


xt b t n ui satisfies (EH| for 

□ 


4 Runge approximation property 


In this section, we derive the Runge approximation property for the following 
anisotropic Maxwell equation 

( V x E — ikpH = 0 

< m £2, 

[VxH + ikeE = 0 

where p, is a smooth scalar function defined on £2 and e is a 3 x 3 smooth positive 
definite matrix. Recall that 


p{x) > po > 0 and E (*)£<& > £ ok | 2 VC € 

*-j=i 


If we set u = 


then we have 


H 

E 


and 


L := i 


c -i 


0 Vx 

0 p~ l h J V ~^x 0 


Lu = 0, 


fc/fi 


where Ij means j x j identity matrix for j = 3, 6 . 


(4.1) 

(4.2) 


Theorem 4.1. Let D and £2 be two open bounded domains with C°° boundary 
in R 3 with D <s £2. If u £ ( H(curl,D )) 2 satisfies 


Lu = 0 in D. 


Given any compact subset K C D and any e > 0, there exists U £ ( H(curl , £2 )) 2 
such that 

LU = 0 in £2, 

and \\U - u\\ H (curl,K) < e, where \\f\\ H (curi,n) = (ll/IU 2 (n) + \\curlf\\ L 2 ( Q) ). 

Proof. The proof is standard and it is based on weak unique continuation prop¬ 
erty for the anisotropic Maxwell system L in 63 and the Hahn-Banach theo¬ 
rem. The unique continuation property of the system L is proved in HU- For 
more details, how to derive the Runge approximation property from the weak 
unique continuation, we refer readers to m- □ 
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5 Proof of Theorem 1.1 


In this section, we want to use the Runge approximation property and the OD 
solutions to prove Theorem 11.11 We define B to be an open ball in R 3 such 
that lie B. Assume that C R 3 is an open Lipschitz domain with B C ft. 
Recall we have set ui £ S 2 and { 77 , w} forms an orthonormal basis of R 3 and 
tg = inf^gu x ■ to = Xq ■ to, where xq = xq (w) £ dD. 


5.1 Penetrable Case 

For the anisotropic Maxwell’s equation 

V x E = ik[iH 
V x H = —ikeE 
div(e£) = 0 
di v(/iF) = 0 , 


(5.1) 


for any t < to and 77 > 0 small enough, in section 3, we have constructed 


Et-rj — 

Ht-r, = 


F\(x)e iTX '^e~ T ^ x ' ul ~ ( ' t ~ 71 ^ A ^ ( ' x, ' > b + 
F\(x)e iTX '^e~ T ^ x ' ul ~ < ' t ~ v ^ A * ( - x '' > b + 


T i£b,t-r l ,N, U ( x > T )+ r ££b,t- 

rT ) + r t 2 b,t- 


r/,N,uj 

rj,N,uj 


{x,t) 

(x,t) 


to be the oscillating-decaying solutions satisfying EH in Bt- v {to) = BC] {x\x ■ 
u > t - 77 }, where F\(x) = 0{t) and F\{x) = 0(t 2 ). Moreover, 
and r xub,t- v ,N, u satisfy EH for M = 1 and |a| = 2, respectively, r^ A b t _ v N bJ 
and r AA b t N u, satisfy (13.81) for k = 1 and k = 2, respectively. Similarly, we 
have 


E t = F\(x)e iTX 'Ze- T ( x ' u ’- t '> A t( x ''>b + 

H t = Fl(x)e iTX 'Ze- T( ~ x ' ul - t '> A ? ( ~ x ''>b + 


^Xt,b,t,N,w( X ’ T ) + r X t\b,t,N,uj(. X ’ r )> 
T tAt,N,J X ’ T ) + r tub, t,N,<A X > T )> 


so be the oscillating-decaying solutions satisfying EH in B t (oj) = B^{x\x-uj > 
t}, where T AA btNuj and T A f b t Nu] satisfy EH for M = 1 and |a| = 2 , respec¬ 
tively, x AA b t N u and x^ A b t N u satisfy (13.81) for k = 1 and k = 2, respectively. In 
fact, from the construction the oscillating-decaying solutions and the property 
of continuous dependence on parameters in ordinary differential equations in 
section 3, it is not hard to see that for any r, 


Et-r] —> E t 
Ht-r] Ht 


in H 2 (B t (co)) as 77 tends to 0. 

Note that fl t (w) C B t -n(uj) for all t < i 0 . By using the Runge approximation 
property, we can see that there exists a sequence of functions (E v ^, i = 

1 , 2 , • • •, such that 


En,e —> E t - V 

Hr],i ^ Hf—r] 


in H(curl, B t (uj)), 
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as £ —> oo, where [E n ^, satisfy (15.11) in £2 for all rj > Q,£ G N. Recall that 
the indicator function I p (r,t) was defined by the formula: 

Ip(r,t) := lim lim I e /(T,t), 

rj—>0 t—too ^ 


where 


I v /{r,t) := ikr f (u x H^i) ■ ((A D - A 0 )fy x H Vi t) x v)dS. 

Jan 

We prove the Theorem 1.1 for the penetrable obstacle case. For the anisotropic 
penetrable obstacle problem 

V x E — ik^H = 0 in fl, 

< V x H + ikeE = 0 in fl, (5.2) 

v x H = f on d£l, 

where k is not an eigenvalue of (15.21) . Moreover, we assume /i is a positive smooth 
scalar function, e = eofy) — xd^d{x), where 70 is symmetric positive definite 
smooth matrix, £d(x) is a symmetric smooth matrix with dete^fy) fy 0 \/x G D 

{ 1 x £ D 

. Moreover, we need e = e(x) is a positive definite 

0 otherwise 

matrix satisfying the uniform elliptic condition. Recall that when e(x ) = £ 0(^)5 
we have constructed E t and H t which are oscillating-decaying solutions defined 
on the half space for the anisotropic Maxwell’s equation 


V x E — ik^iH = 0 in fl, 

V x H + ikeE = 0 in il. 


(5.3) 


and {(E Vi t, H v ^)} are sequence of functions satisfying (15.31) defined on the whole 
Q. Therefore, we can define the boundary data f n j. = v x H rj j on dfl and solve 
{E,H) satisfies (15.21) . Let H n j = H — H v 7 be the reflected solution, then 
satisfies 

JV x (e -1 V x 7 ) — k 2 ^H v ^ = —V x ((e - 1 (:r) — £q 1 (a;))V x H^f) in f l, 

1 v x H r] j = 0 on 

(5.4) 


Lemma 5.1. We have the following estimates 
1 . 

> f [e(e -1 - £o 1 ) _1 £o ' 1 V x H V<1 \ ■ (V x H v ^)dx - k 2 f [i\H V: e\ 2 dx. 
J d Jn 

2. 


r 1 I*' £ (T,t)> / ((e^-e 1 )V x H Vf t) • (V x H^ijdx - /c 2 / n\H lh z\ 2 dx. 


ID 


Proof. First, we need to prove the following identity 

-r-fy^fyfy) = [ ((e - 1 - efy 1 )V x H Vtt ) • (V x H^)dx 

Jn 

- f (e -1 V x H~!j) ■ (V x H~i,i)dx - k 2 f /i|ify^| 2 ((5.5) 
J D, J D 
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Multiplying H v i in the equation (15.41) and integrating by parts we have 

! [ 

in 


(e i V x H Vt i) ■ (V x H Vi t)dx — k 2 


i dx 


/(( 


+ / ((e " 1 - £q')V x 77^) • (V x H^)dx = 0, 


in 


/ (e _1 V x H v/ ) ■ (V x H v/ )dx - k 2 
in 

- / (e _1 - e^ 1 )V x 77^) • (V x H Vtt )dx 
in 

= - [ ((e -1 - £o X 77.^) • (V x TJ)dx. 
in 

On the other hand, H{x) satisfies 

V x (e _1 (x)V x H(x)) - k 2 fiH(x) = 0, 


dx 


(5.6) 

(5.7) 

(5.8) 


then multiply by 77 TJj ;( x) in the equation (15.81) and integrating by parts we have 


/<«-■ 

in 


- ey 1 )V x H rh t;) ■ (V x H)dx = / (e _1 V x 77) ■ (y x T7„,*)ds 


'an 


/ 

Jan 


V x 77,.,,*) • (i/ x T70&0) 


Thus, combine (15.61) . (15.91) and f dn (v x 77,,,*) • (e 0 1 V x H v ^)ds is real, then we 
have 


/<«-■ 


V x 77,,,*) ■ (V x H v ^)dx — k 


I 

in 


dx 


- [ (( e_1 - £o : )V x if,,/) • (V x 77„,*)dz (5.10) 

in 

[ {vx 77,,,*) ■ (e _1 V x 77)ds - [ (y x 77) • (e^ V x 77,,,*)ds 
JdQ Jdfi, 

- / iy x 77 r ,,*) • (e _1 V x 77)<is — / (i/ x 77,,,*) • (cq 1 V x H v ^)ds 

Jan Jan 

- / iy x 77,,,*) ■ (e _1 V x 77)ds — / (i/ x 77,,,*) • (cq 1 V x H v ^)ds 

J dfi J dfl 

- / (^ x 77,,,*) ■ [—ikE + 7fc75,,,*]ds 
Jan 

-ik / (u x 77,,,*) ■ [(Ad - A@)(i/ x 77,,,*) x ^]ds 
Jan 


=T~ 

P 

Second, we show the following identity 


/ (^o 1 

in 


V x 77,,,*) ■ (V x H^^dx — k 2 / /r|i7^^| 2 da; 


(5.11) 

(5.12) 


+ / ((e _i 0r) - £q i (x))V x 77) • (V x i7)dx 
in 

= -r~ l r/. 
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Replacing H r] i(x) by H(x) — H rl j(x) in the equation (15.41) . then we have 
V x ((e ” 1 - Co 1 )V x H) + V x (e^V x H~~ f ) - k 2 ^H^e = 0 in fi. (5.13) 


Multiplying H v j(x) in the equation (15.131) and using integration by parts we 
have 


j f ((e " 1 - eu 1 )V x H) ■ (v x 


dx 


+ / (e 0 V x H rhi ■ (V x H r]t \dx-k fi 


) dx - k2 L 


ii, 


ij,£ 


dx = 0 , 


since v x H v i = 0 on <917. Then we can write equation (15.141) to be 

J 1 V x H rlJ ^j ■ x dx — k 2 fj, 

+ [ ((e" 1 - Cq X )V x H) • (V x H)dx 

J n 

= [ (( e_1 - e^)V x H) • (V x H^)dx. 

Jn 

Eliminating H(x) by H v j(x) + H v j(x) in (15.151) we have 

J (cq x H n j^j ■ (V x dx - k 2 J n 

+ [ ((e" 1 - 1 )V x H) • (V x H)dx 

Jn 

= [ ((e -1 (x) - Cq 1 (x))V x H v> e) ■ (V x H^^dx 

Jn 

+ [ ((e _1 (s) - e ( )' 1 (a;))V x H Vt i) • (V x H r] ^)dx 

Jn 

Again from (15.41) and by taking the complex conjugate, we can write 


(5.14) 


(5.15) 


(5.16) 


V x (e -1 V x Hr, t e) - k 2 ^H v/ + V x ((e~\x) - £q i (x))V x H^) = 0. (5.17) 

Multiplying by H V: i(x) in the equation (15.171) and using integration by parts we 
have 


(e X V x H Vt t) ■ (V x H Vt e)dx — fc 2 / fj,\H Vi i\ 2 dx 


+ 


f ((e-'C 

Jn 


x) — e 0 i (a;))V x H^j) • (V x H v ^)dx = 0. 


(5.18) 
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Then from the equations (15.161) . (|5.18l) and the first identity (15.51) . we can obtain 
/ (fig 1 V x H Vtt ) • (V x H Vti )dx - fc 2 / n\H v .e\ 2 dx 

+ f ((e _ 1 (:e) — e^" 1 (x))V x H) ■ (V x H)dx 
Jn 

= [ ((e - 1 (:r) - e ( 7 1 (x))V x H Vi t) ■ (V x H Vt e)dx 


j 

- / (e _1 V x H v>i ) ■ (V x H Vt i)dx + k 2 n\H v> i\ : 

J r2 Jo, 


dx 


in 

= - 


(5.19) 


Combine (15.191) with the formula 

(cq X V x H^) • (V x H^) + ((e” 1 - Eg : 1 )V x H) ■ (V x 7?) 

=((e _1 - Eg : 1 )V x H) • V x H + Eg 1 (V x H) ■ (V x H) 

- 2Re {eq 1 V x H • V x } + Eg 3 V x H r] ^ ■ V x H r] j 
=e _1 (V xif).(Vxfl)- 2Re {e^V xH-Vx 77 ~i} + e^V x H fl ■ V x 77^ 


?V x if- e^e " 1 (V x H^ e ) ■ e~? V x H - 1 (V x H v/ ) 

eh^ 1 (V x 77^7) 


£ 2 e^ 1 (V x H. nJ ) 


Eg Vx H^£ ' ^ X 


■5V x H-eho 1 (V x H^e) ■ e~iV x H - eh^ 1 (V x H v ^) 

+ (fig 1 — ee o 2 ) x • (V x 77)77) 

>[(/ - ££g 3 ) Eq'V x H v/ ] • (V x 77)77) 

>[e(e _1 - eo 1 V x H^\ ■ (V x 77)77) 

and note that 


e'sVxff- e^e" 1 (V x 77)77) • e~* V x i? - e^e” 1 (V x 77)77) 

Therefore, we get 

[ [e(e -1 - eo" 1 )- 1 '- 1 - - - 1 ^ J “ " 2 

J D 


> 0 . 


-T^lf > 


L V x H v> t] ■ (V x H Vit )dx - k“ 


[ lAHr],l\ 

Jn 


dx 


which finished the part 1 of lemma 4.1. Finally, again from (15.111) . we have 


> 


/((. 


1 - t 3 )V x H rh t) ■ (V x H Vi i)dx — k 2 I y\H v ^{ 2 dx. 


□ 


Remark 5.2. The first inequality will be used when (e 1 — e 0 3 ) is strictly pos¬ 
itive definite, i.e. 

£ • (e ” 1 — Eq *)£ > A |£| 2 for all (el 3 and for some A > 0; 

and the second inequality will be used when (eq 1 — e - 1 ) is strictly positive 
definite, i.e. 

£ • (e ^ 1 — e _1 )£ > A |£| 2 for all (el 3 and for some A > 0. 

Now, our work is to estimate the lower order term H v ^. 
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5.1.1 Estimate of the lower order term H v ^ 

Proposition 5.3. Assume fl is a smooth domain and D <s 12. Then there exist 
a positive constant C and 6 > 0 such that 

\\Hr,,e\\L 2 (Q) < C||V x H V} e\\ L p(D) 

t t 2 + , 

for every p £ (max{ -, —— }, 2 . 

3 l + o 

Proof. We follow the proof of the proposition 3.2 in [?]. Fix l £ N and we set 
/ := ~(e -1 - to 1 )(V x H^j), g = 0. Note that, e -1 - e+ 1 = e~ 1 (e D XD)eo 1 is 
supported in D. Then the reflected solution H v ^ satisfies 

Jv x (e _1 VF^) — k 2 fiH Vt £ = —V x ((e _1 (j;) — e ( ^ 1 (x))V x in 12, 

1 v x H n j = 0 on 512. 

(5.20) 

From the L p estimate (Theorem 6.6), if we consider the following problem 

f V x (e^V x U) + e-\JJ = V x / in 12, 

| v x U = 0 on dll, 


has a unique solution in Hg' q (curl, 12), where e^ax is th e maximum value among 
all eigenvalues of the matrix e _ 1 (x) in the region 12. Moreover, we have the 
estimate 


\\U\\ LP{ n) + ||V x U |Up ( o) < C\\f\\ LP{p) (5-21) 


for p £ (y——, 2 ] for some S > 0 which depends only on 12 . 
II = H Vt e — {7, then II,^ satisfies 


Now, we set 


I V X (e 1 vn^e) - k 2 pn v>£ = (k 2 p + eJ ax )U in 12, 
x n^ e = 0 on <912. 

By the well-posedness of (15.2211 in H(curl , 12) for the anisotropic Maxwell’s equa¬ 
tion (see Appendix), we have 


l|n»)/llL 2 (n) + ||V x n^^|| L 2 (n) < C'||{ 7 || i 2 ( fJ ) 
if k is not an eigenvalue. Moreover, for p < 2, it is to see that 


(5.23) 


]|n^^||x,p(n) + ||V x n^y|| iP (Q) < C\\U\\ L 2 (ny 
Following the proof in the proposition 3.2 in [7] again, we denote B p £ 2 (Vl) to 

P 

be the Sobolev-Besov space, then we have U £ B p f (12) and the inclusion map 

V 

B p £ 2 { 12) —» L 2 { 12) is continuous for p £ (|,2]. Moreover, since V x U = 0 and 

P 

v x U = 0 on <912 and use Lemma 7.6 ( property 5 in the appendix of Cl) , we 
have the estimate 


\\U\\l*( n) < CWII B p.’ (n) < C{\\U\\ LP(a) + ||V x 27||z,p ( n)} (5.24) 

P 
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for p £ (|, 2]. Combining (15.211) . (15.231) and (15.241) . we obtain 

||n^|| L p ( n) + ||V x II J) . f || i p( n ) < C'||/|| i p(n) (5.25) 

for p £ (max{|, |±|}, 2]. Since H v j = 11,,^ + U, by using (15.211) and (15.251) . we 
have 

I|£^IUp(Q) + ||V x H^\\ LP{n) < C||/||iP (n ). (5.26) 

Since v x H v i = 0 on dSl, we use the Lemma 7.6 again, then we can obtain 

\\H v ,i\\L 2 (n) < C||^r),^|lsp. 2 (n) 

P 

< C{\\H Vt t\\ LP (ty + ||V x H Vt ( IIlp(q) + || V ■ iL I?i f|| i p(n)(|i.27) 

In addition, from (15.201) . it is easy to see 0 = S7-(p.H Vt e) = V/j ■ H v ^+p(V ■ H^j), 
then we have 

||V • ^IUp(O) < "7^ 11^11^. (5-28) 

Finally, use (I5.26|) . (15.271) and (15.281) . we will get 

11 # 77 , f|U 2 (n) < C'{||J# TJ ,^||z,P(n) + || v x #^|| L p(fj)} 

< c'II/IUp^) 

< C\\X7 x H v , e \\ LP(D) . (5.29) 

□ 


Remark 5.4. In the reconstruction scheme, we need to take limsup^go for 
(15.291) on both sides and i7 t+?? —>• H t in H{curl , Sl t {oj)) as p — > 0, then we have 


lim limsup 11^^11^2(0) < C||V x H t \\ LP i D ), 


for p £ ( — ,2]. Moreover, if (15.201) is written as the following form 
3 


V x (e^VH) - k 2 pH = -V x ((e" 1 ^) - 1 {x))V x H 0 ) + k 2 (e - e 0 )H 0 . in SI, 

v x H = 0 on dSl, 


and we can derive the following estimate by using the same method in the proof 
of the Proposition 5.3, then the estimate (15.291) will be 


||#||L 2 (n) < C{||V x #o||lp(d) + ||#o||L 2 (n)}, 

for p £ (^, 2 ], 

In view of the lower bound, we need to introduce the sets Dj } $ C D, Dg C D 
in the following. Recall that hp(p) = ■ p and to = h]j(p) = xo ■ p 

for some xo £ dD. Va £ dD D {x • p = 1 id{p)} '■= K, define B{a,8) = {x £ 
K 3 ; \x — a\ < 5} (S > 0). Note K C U Q ,£kB(cx,S) and K is compact, so there 
exists ai, ■ • • , a m £ K such that K C UjL 1 B(aj,S). Thus, we define 

Djj := D D B{aj,8) and Ds := U JL 1 Dj > g. 
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It is easy to see that 

jI d\d s e- pr{x -“- to) <^bdx = 0(e~P aT ) 

1 / d \ z ? 5 e- pTix - u - to) < ix,) bdx = 0(e~P aT ) 

where Af- Q ( x '), Af ( x ') are smooth matrix-valued functions with bounded entries 
and their real part strictly greater than 0. so 3a > 0 such that ReA^ (x') > a > 0 
and ReA^(x') > a > 0. Let ay- S K 1 by rotation and translation, we may 
assume ay = 0 and the vector aj — Xq = —Xo is parallel to e3 = (0,0,1). 
Therefore, we consider the change of coordinates near each aj as follows: 

fy' = x' 

\z/3 =X- p-to, 

where x = (xi,x 2 ,x 3 ) = (x',x 3 ) and y = ( 2 / 1 , 2 / 2 , 273 ) = ( 2 /', 2 / 3 )- Denote the 
parametrization of dD near aj by lj{y'), then we have the following estimates. 
Note that the oscillating-decaying solutions are well-defined in D. 

Lemma 5.5. For q < 2, t 1, we have the following estimates. 

1 . 


\H t (x)\ q dx < 


id 


T 2q ~ l Yl // e- aqTl ^ v,) dy' + 0 (T 2 q - 1 e~ qa5T ) 

]ri-lj\y'\<t> 

+0(T 2q e~ qaT ) + 0(re~ CT ) + 0(t~ 2N+5 ) 


2 . 



> 


Ct 3 // e~ 2 aTl ^ y,) dy' - Cr 3 e~ 2aST 
-CTe~ 2cT - Ct ~ 2N+5 


3. 


h 


\E t (x)\ q dx < 


m .... 

r^V // e- aqTl ^ v,) dy' + 0{T q - 1 e~ qaST ) 

j =1 JJ\y'\<6 

+ 0 (r q e~ qaT ) + O (r- 1 ) + 0 (t~ 2N+3 ) 


4■ 



> 


m .... 

Cry'// e- 2 aTl ^ y,) dy'~CTe~ 2aST 
J=1 JJ\y'\<S 

-Ct” 1 - Ct ~ 2N+3 


Proof. The proof is via the representation of the oscillating-decaying solutions 
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of (E t ,H t ). For r l(r <C r 2 ), we have 


f \H t \ q dx < Cr 2q f e- qaT ^~^dx + C q f \T^ N J q dx 

JD JD JD 


+c ( 


L 


A, 2 

r Xt,b,t,N,u 


| q dx 


< Cr 2q [ e - qaT (z-u-*o)dx + CT 2q [ e- qaT{x - u ~ to) dx 
Jd 5 J D\Ds 

+ C q [ l r AB^^dx + Cq [ \r l A ,B, 1 J q d x 
J D J D 


m f'f* r^ 

< Cr 2 q y\ // dy' e- qaTy 3 dy 3 + CT 2 q e~ qaT 

j =l JJ\v'\<& 

+ C\\r^ 2 btNt J\ 2 L 2^D) + C \\ r x^b,t,N,u\\ 2 L 2 (D) 

< Ct 2 ^ 1 V [[ e - aqT l^y') d y' - — T 2q - 1 e -qaSr 

~iJJ\y '\<6 9 

+CT 2 q e~ qaT + Cre- caT + Ct~ 2N+5 , 

where c is a positive constant and a depending only on a A ,aB- For the lower 
bound of f D \H t \ 2 dx, we have 


j D \Ht\ 2 dx > CT4 ^ e “ 2aT(K '“ _to) ^- C ll r x; 2 M,iv,Jli=(n f0 


ID 

A,2 


W) 


C \\ r Xt,b,t,N,J^ 2 {n t0 (u)) 


> Ct ] 
Ct 3 


f 

JDs 


- 1 aT{x-u-t 0 ) dx _ c T e~ CT - Ct~ 2N+5 . 


lib rtrt 

V // e~ 2 aTl ^dy' - CT 3 e~ 2aST 

~[JJ\y'\<b 

-Cre- caT - Ct~ 2N+5 . 


It is similar to prove the remaining case, so we omit the proof. 
Lemma 5.6. IFe have the following estimate 

m\i> {D] 


□ 


ii^ 


p- > 0(t ), r > 1. 

11 LHD) 


Proof. Since dD is Lipschitz, we have lj(y') < <I7|j/'|. Therefore we have the 
following estimate 

Ct 3 jr [[ e- 2aH MU y ' > Ct 3 ]T // e~ 2a ^'\ 


o =i 


=i JJ \y'\<b 


3=1 


=i JJ \y'\<b 


> 


Cr E // e- 2 ^V 

■j—l JJ\v'\<tS 

0(t). 
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Then we use lemma 4.4 to get 


\mh (D) 

ll^ll| 2 (D) 


1 - 


Ce~ 


r +Cr- 


t +Ct~ 


Ct 2 


If\y'\<6 ' 




dy' 


0(e~ 2l5a ' r )+0(Te~ caT )+0(r ~ 2Ar+2 ) 
T.f=i ffiy'Ks e~ 2a ’ rl ^ y ' ) dy' 


1 - 


= 0(r 2 ) (if r > 1). 


□ 


Lemma 5.7. If t = Hd(p), then for some positive constant C, we have 

liminf / r|V x H t \ 2 dx > C. 

T ^°° Jd 

Proof. Since lj(y') < C\y'\, we have 

I |V x H t (x)\ 2 dx > C I \E t (x)\ 2 dx 
Jd J d 


Therefore, we have 


> CrV // e- 2 aTl ^ y,) dy'-CTe~ 2a6T 

J = 1 JJ\y'\<6 

-Ct- 1 - Ct~ 2N+3 

m pp 

> e- 2aT ^Uy' ~CTe~ 2aST 

]=1 JJ\y'\<5 

-Ct- 1 - Ct~ 2N+3 

m pp 

> Ct[t~ 2 ^2 e- 2a \ y '\dy'}-CTe- 2aST 

j -1 JJ\v'\<tS 

-Ct- 1 - Ct~ 2N+3 (as r » 1). 


liminf / r|V 

T ^-°° Jd 


x H t \ 2 dx > C. 


□ 


Lemma 5.8. For p £ (max{|, y±|},2]. we have the following 


lim lim sup 

i->oo 


\\V * Ht\\ 2 L2{D) 


< Ct 1 - 


2 

P 


(t » 1 ). 


Proof. From the proposition 5.2, we have 

lim limsup 11^^11^2(0) < CUV x H t \\ LP{D) . 

r/->0 


Then it is easy to see the conclusion. 


□ 


Remark 5.9. Recall that the sequence {H v e} converges to H t+V in H(curl, K) 
as l —> oo for all compact subset D <s K <s 14 and H t + V —> H t in if 2 (fi t (u;)) as 
77 —> 0 , so we have 

||V x H v j ||lp(_d) -t ||V x H t ||lp(d) and \\H v ^i?( D ) ->■ \\H t \\ L 2 ( D ) 
as I —>■ 00 , r) —> 0. 
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5.1.2 End of the proof of Theorem 1.1 for the penetrable case 

First, we prove the case t < ho{p). From (15.51) . we have 





6 q )V X • (V X H^^dx 

V x Hrj^) ■ (V x H v ^)dx — fc 2 / p\H v< (\\d>x}>Q) 

in 


Note that ( E e H t j) satisfies 


V x E r i,f. — ikpH v ^ = 0 in Q, 

V x H Vt g + = ik(e o — E)E v j in fl, 


and rewrite it as 


V x (e x E v>e ) - fc 7 -E^ = k 2 (e - e 0 )E 7l/ . 


(5.31) 


Thus, we can use the same argument from the Remark 5.4 again to (15.3111 . it is 
easy to see 

\\E v ,e\\L 2 (n) < C\\E v ,t.\\ L 2( D )- 

In addition, we use the Maxwells equation and e — eg = —cdXd, then we have 




/( 


V x • (V x Ed r) ^)dx = / (— ikeE v ^ + ik(e o — e)E r/t e)) • (V x H^^dx 


C 


f \Erj,e\ 2 dx + C f 

in i l 


dx + C I \E v/ \ 2 dx (5.32) 

D 


L 


C I lErff^dx. 


Thus, from (15.301) . Proposition 5.3, Lemma 5.5 and (15.321) . we can obtain 

\~kp’ e (T,t)\ < + \\Hr],e\\ 2 H(curl,D)- 

From taking t —> oo and 77 —>■ 0, we have 

l lp[T,t) | < |T ^ / 

j= 1 J J \v'\<& 


m .... 

,(r, t)| < | tV e- 2aTl ^ v,) dy' + 0{T 2 e~ 2aST ) 


In particular, we get 


+0(r 2 e~ 2aT ) + 0(r~ 3 ) + 0 (t~ 2N+3 ) 
< 0(t _1 ) +0(r 2 e~ 2aST ) 

+0(r 2 e~ 2aT ) + 0(r~ 3 ) + 0 (t~ 2N+3 ) 


lirnsup |—/ p (t, f)| = 0. 

T—> 00 T 


Second, we prove the case t = Hd(p)- 
Case 1. £ • ( 7” 1 — 7 c j“ 1 )£ -2 A |£| 2 for all (€l 3 for some A > 0. 
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From the inequality in Lemma 5.1, we have 


1 TV,* 
P 


> 


h 


) 1 e 0 : V x H Vi i] ■ (V x H Vi i)dx — k 


/ k\ 

J a 


2 dx 


> 


-k 2 [ nlH^^dx 
■In 

C JjV x H^fdx - c\\H^e\\ 2 LHny 


By using the definition I p (r,t) := lim^o lim<-too /p’ f (r ,t), {H v /\ converges to 
H t in H(curl, K) for all compact subset D <s K (s ft as £ —> oo, T] —> 0, we have 


~ I p( T , t ) 

l|Vxif t ||’ a(D) 


Ct 


1 — C lim lim sup ■ 

£_>0 ^->oo 


™ 2 


L 2 (n) 


IV x H t 


'L 2 (D) 




Hence, using Lemma 4.7 we deduce that for t ^ 1, 


I I p {t, h D {p)) | > C > 0 


which finishes the proof. 

Case 2. £ • (y ^ 1 — 7 _1 )£ > A |£| 2 for all £ G M 3 for some A > 0. 

Similarly, using the inequality in Lemma 4.1, we have 

T~ l P p ' l {r,t) > f ((e^ 1 - e _1 )V x H V} e) ■ (V x H Pt f)dx - k 2 [ p\H v ^\ 2 dx. 
J d J n 

Then use the same argument as in Case 1 we can finish the proof. 


5.2 Impenetrable Case 

We give the proof of the second part of Theorem 1.1, since it is the hardest 
part. The other cases are easy since we have proved it in the penetrable case. 
In addition, the upper bound is easy because of the wcll-posedness and the L p 
estimate for the indicator function, but the lower bound is not easy to see. In 
the following proof, we will use the layer potential properties for the exterior 
isotropic Maxwell’s equation (with the Silver-Miiller radiation condition) and 
the perturbation argument from the anisotropic Maxwell’s equation compared 
with the isotropic case. In the impenetrable case, we have chosen the oscillating- 
decaying solution as the following form 


Et = G 2 B (x)e iTX <e- T ^ u - t ^^b + T^ 2 btNt 
H t = G] } {x)e iTX <e- T ^- t ^ A ?^^b + T^ tN 


( \ . ,B, 2 

( X >T) + r xt,b,t,N,u 

(X, r) + r^ b>tiNtU 




where G B (x) = 0(t) and G 2 B (x) — 0(t 2 ) and V^f bttNul satisfies (Id.91) for 
|a| = j and r^’ J b t N u satisfies (13.91) for k = j. 

We start by the following lemma. 

Lemma 5.10. Assume that fi is a smooth scalar function and 7 is a matrix¬ 
valued function. Let (E,H) G H(curl; Q\D) x H(curl; fl\D) be a solution of 
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the problem 


' V x E — ikpH = 0 in fl\D, 

V x H + ieE = 0 in fl\D, 

v x E = f on dTl, 

v x H = 0 on dD , 


(5.33) 


with f 6 TH 1 / 2 (9J7). If we put f Vl e = v x E v j with {E v j} is obtained by the 
Runge approximation property. Then we have the identity 




- [ (|V x E Vtt {x )| 2 - fc 2 |^T,,f(ar)| 2 }da; 
Jd 

- f {|V x £^(a;)| 2 - k 2 \E^e(x)\ 2 }dx 
Jq\d 

f {|V x H Vt i(x)\ 2 — k 2 \H 71 ^(x)\ 2 }dx 
J D 

+ [ {|V x H v ^(x)\ 2 — k 2 \H V} t(x)\ 2 }dx 

Jn\D 


and the inequality 


1 


—iy{r,t) > / {|V x HM\ Z - k 2 \H^(x)\ 2 }dx - fe- 


' D 


/ j 

Jn\D 


\H v ,e(x)\ 2 }dx, 


where E ^ = E — E^^ and H v £ = H — are described in section 5. 

Proof. Use the integration by parts and the boundary condition, we have 

I e _1 (Vx£HV x E^e)-k 2 eE-E^ e dx =-( [ -[ )ik{uxH)-E^ e dS = 0. 
JQ\D Jd f2 J dD 

Adding this to 

Ip’^ = [ (v x E v ^) ■ (—ikH + ikH V 'i)dS 

Jd n 

= / — (m _ 1 V x E Vt t) ■ (V x £) + k 2 {p,E v ^) ■ Edx 

Jn\D 

+ / /x _1 |V x -E^l 2 — k 2 {pE v ^) ■ E Vt edx + [y x E Vt e) ■ (-ikH)dS 
J O ’ 

due to the zero boundary condition on dll we have the last term is vanishing. □ 

Fr om the above estimate, it only need to control the lower order term 

fn\B \ < J jX - 


5.2.1 Estimate of the lower order term H v ^ 

Proposition 5.11. Let fl be a C 1 domain, Dj§ Lipschitz. Then there 
exists a positive constant C independent of (E Vt e, H Vt i) and (E v> e,H Vt i) such 
that 

[ \Eri,e{x)\ dx < C{||V x H V: (\\ LP ^ D ^ + \\H V! e\\ Ha+1 / 2 , D s}, 

Jn\D 

for all p and s such that max{2 — <5,4/3} < p < 2 and 0 < s < 1 with S > 0. 
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Proof. Step 1 . Before proving the Proposition 6.2, we consider the anisotropic 
Maxwell’s equation in 12 as follows: 


! V x E r] j — = 0 in 12, 

V x H Vi £ + ikeE v ^ = 0 in 12, 

v x E rlt f := f rj i £ Tff _1 / 2 (912) on <912, 


(5.34) 


E, 


T),t 


are solutions of the anisotropic Maxwell’s equation. Since 
> = H — we have 

f V x E V} e — ikjiHr= 0 

in 12 \D, 

1 V x H Vt e + ik^Er-jj = 0 

“ »n 5 ’ < 5 - 35 ) 

| v x E Vt i = 0 

on 912, 

Z/ X £ — V X l 

on dD. 


Step 2. Let {E^,E[^ e ) be the solution of the following well posed exterior 
Maxwell’s problem 


'V x E% - ikH% = 0 
V x + ikE™ = 0 
v x H** e = -v x H Vtt 

satisfiy the Silver-Miiller radiation condition. 


in R 3 \Z2, 


in M. 3 \D, 
on dD 


We can represent these solutions E and H^ x t by the following layer potentials 
H^(x) := Vx/ ®k(x,y)f{y)ds(y), 

J dD 

E™Ax) := --J-V x m x A x), x £ R 3 \dD, 

h ik " 

e ik\x-y\ 

where &k(x,y) = — -—:-r, x,y £ R 3 , x + y, is the fundamental solution of 

47r|a; — y\ 

the Helmholtz equation and / is the density. Now, we follow the arguments in 
section 2.1 of [7] and use the same argument for the isotropic Maxwell’s equation 
(15.361) . then we have 


\\Ef 1 fi\\LP(a.\D) < C{\\v x H Vt e\\ L P(dD) + ||V x ||lp(d)}, ,g 

ll^/llL 2 (n\S) < C{\\u x H r] ^\\ L p(dD) + ||V x H Vt e\\ L p( D )}, 

for p £ (-,2]. Moreover, if we define £ rh t; = E v ^ - E^ e , 'H v ,e = H rj j — Hf-f , 
then £ V} g and satisfy the following Maxwell’s equation 

{ V x £ v>t - ikpH, h e = ik( 1 - 
V x Hr,,e + ike£ Vt j( = ik(j - h)E^ t 
v x 'H v ,i = 0 
v x £ Vti = -v x E^ t 


in 12\Z), 
in 12\Z), 
on 912, 
on dD. 


(5.38) 
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Step 3. Now we decompose £ v j, = £* t + £ 2 t and 77 r ,y = 77^ + 77 2 ^, where 
(£^ti'Hlit) satisfies the following zero boundary Maxwell’s equation 


f V x = ik( 1 - n)H% in Q\D, 

< V x + ikt£\ t = ik(e - h)E “ in 12\A 

x = 1 / x 77^ = 0 on d(£l\D), 

and (£ 2 ^,77 2 J satisfies 


{ V x £^ e - ikpH 2 ^ = 0 in f2\.D, 
V x 77 2 ^ + ikj£ 2 e = 0 in £l\D, 
v x ^r,y = 0 on d£l, 

v x £^ e = -v x E^ J j on 377. 


(5.39) 


(5.40) 


First, we deal with the equation (15.3911 by using the L p estimate in tt\D. 
Note that (£* t ,7 7* ^) satisfies (15.391) . then we have 


V x (e 'Vx £^ e ) - k 2 = ikV x [(/i 1 - 1)71®^] + ik (7 - h)E ^ 

^ x = 0 


in fi\77, 
on d(£l\D), 


and 


V x (e " 1 V x 77^) - fc 2 /i77^ = ik\7 x [(/ 3 - e" 1 )^] + ik{ 1 - /x)77^ 
^ x ^ = 0 


in f2\77, 
on d(n\D). 


Now, if we use the same method in the proof of the Proposition 5.3, we will 
obtain 


| W^,e\\LP(n\D) + ll ^ 7 x ^,e\\LP(n\D) < ^{\\^^A\lp(q\d) + ll-^e/lli/qnyD)}) 

\ll^^llip(n\B) + l|V x ^i,i\ U p (n\B) < C{ll- E '^llLp(n\B) + ll^/lli 2 (n\5)- 

(5.41) 

4 

for any - < p < 2. If we combine (15.3711 and (15.4111 together, we have 
o 

llip(n\B) < C'dl*' x Hrf^W L p(OD) + IIV x H v ^\\ L p^ D )}. (5.42) 

For (£ 2 ^,77 2 (), we apply the L 2 -theory for the anisotropic Maxwell’s equa¬ 
tion, we get 

llz, 2 (n\D) < \\£ri,A\H(curl,n\D) < C\\vx£ 2 (\\ H -i/ 2 ( a n ) < C\\l/xE^ X e \\ H ^i/ 2 ^ d Q) . 
Moreover, following the proof in the Lemma 2.3 of [7j, we have 


\\ v x E^ e \\ H -i/2(dn) < C\\f\\LP(dD), Vp > 1, 

and 


WKA\l^\d) < x ^,/HL.(fl D ) + IIV X H v 4 2 lhd) }, (5.43) 

4 _ _ 

for all p £ (—, 2]. Recall that 7 i v ,e = 77*^ + 77 2 ^, by using (15.4211 and (15.431) . 

O 

then we have 


ll^/IL 2 (n\D) < C{||^ x H Vt e\\ LP ( dD ' ) + ||V x Hrj t e \\lp(d)} (5.44) 
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4 

for all p £ (—, 2]. Finally, for s > 0 and p < 2 we have H s (dD) C L 2 {dD) C 
o 


L p (dD ), then we reduce that 

II v x < C||-ffj)/||LP(a_D) < C||-ff»7^l|j?“(8r>)- 

Note that the trace map from H S+1 / 2 (D) —> H s (dD) is bounded for all 0 < s < 
1. So the estimate (15.4511 will become 



4 

for all p £ (—,2] and 0 < s < 1. 

O 


□ 


Remark 5.12. Now, if we take i —> oo and e —>• 0, we will get 

limlimsup [ \Hr>Xx)\ 2 dx < C{\\H t \\ 2 H . +1/2 , D) + ||V x H t \\ 2 LP(D) }, 

where H t is the oscillating-decaying solution defined on 12 t (w). 

We have the following lemmas for the oscillating-decaying solutions in the 
same way as we did in section 5, so we omit the proofs. 

Lemma 5.13. For 1 < q < oo, r > 1, we have the following estimates. 

1 . 



FO^e-^) + Oir- 1 ) + 0 (t 


-2N+3 


2. 



-Ct- 1 - Ct~ 2N+3 


3. 



+0(T 2q e~ gaT ) + 0(re~ CT ) + 0 (t 


-2N+5 


l 



-Cre- CT - Ct~ 2N+5 
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Lemma 5.14. We have the following estimate 


l-^tlll 2 (D) 


< 0(t 2 ), r > 1 . 


Zt M L 2 (D) 

For p < 2, we have the following estimate 

ll v x H t\\ 2 L P( D ) 

--—j2- < Ct *>, t» 1. 

||V x iJ* || L p^ D ^ 

Lemma 5.15. If t = ho(p),then for some positive constant C , we /lave 


liminf f r|V x U t | 2 <ia; > C. 
r - > °° Jd 


5.2.2 End of the proof of Theorem 1.1 for the impenetrable case 

By using the same argument in the penetrable case, it is easy to see that 

limsup |— I p (t, t )| =0 

t—» oo T 

for i > hjj(p). Recall that from Lemma 6.1, we have 

--If(r,t)> [ {| V x H v ^{x)\ 2 — k 2 \H Vi g(x)\ 2 }dx — fc 2 /" \H^ e (x)\ 2 }dx. 

T Jd Jn\D 

(5.46) 

By using Proposition 6.2, we deduce 

— J^{\^ x Hri,i{x)\ 2 —k 2 \H ri ^(x)\ 2 }dx—C{\\H t \\ 2 HB+1/2 ( D j+\\'V x H t \\ 2 LP ( D j}, 


4 ll^t ll//s + l/ 2 f£)'\ 

where 0 < s < 1 and - < p < 2. We want to estimate —- — .. —-- for 


0 < s < 1. Set r = s + 1/2, then we need to estimate 

\\Ht\\ 2 Hr(D) 


L 2 (D) 


V x H t 


L 2 (D) 


,1 3, 


for r e (-, -]. Using the interpolation inequality, we have 


\H t \\ H r (D ) < C\\H^[ D) \\H t \Y HHD) , 0 < r < 1. 

i + i 

a ' ft 


By the Young’s inequality ab < 5 a — + cF 3 ^-, - + i = 1, we obtain 




C 


5~ a 

a 


Ht\\ L 2 (D) + \ \Ht II // 1 (D) 

{(1 - + r<r- 1 }|| J ff 0 |li 2(D ) + rS r ~' \\VH t \\l{fcfh 


Recall that H t = ^ >, r)+rf; 1 b , tJV Jx, r) 


Xt,b,t,N,u;\ 


is a smooth function with Gg(x) = 0(r) and 1’^’ bt N w satisfies (13.91) for |a| = 1 
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and r^’ 1 bt N u satisfies (13.91) for k = 1. If we can differentiate H t component- 
A dH t dG^e^e-^-^U , dV^ b N dr % 4 

wisely, we will get s ---+ — xt ’”’ + xt ’ 

OXj OXj OXj OXj 

and 

ajut 

| A,B, T.M || 2 < /7_4 f -2 ar(x-p-t)j 

I IIl 2 (d) - or Jd e ax ’ 

BT 2,t 

l^“IU.,D,<cr->A e -. 

< ^-' ,+3/2 . 

Then by using the same method as before, it is easy to see that 

3 o jj 

j =1 J 


< Cr 4 [ e~ 2< ' x ' p ~ t ' > dx + cr _ 1 e _2cT + cr _2Ar+3 . 
J D 


For t = hnip), we have 


II Vff* " 2 


< 

Cr 4 [ e~ 2a{x ' p - hD{p)) dx 

+ ct 4 e 2cT + CT 


J D 




< 

Ct 4 { f + 

[ ) 

g-2 a(x- 

p - h D( p )) dx + CT- 1 


Jd s J 

’d\d 5 




+ct~ 2N + 3 




< 

m pp 

f 

j= i JJ \y 

dy‘ 

1«5 

■f 

e- 2aTV3 dy 3 + Ct 4 , 


+CT~ 1 e~ 2cT 

+ CT~ 

2N+3 


< 

m pp 

L 

j =1 JJ \ V 

e~ 

\<5 

2 arlj (y 

><V - Cr 3 e- 2aST 


+CT 3 e~ 2acT 

+ CT~ 

4 e~ 2cT 

+ ct~ 2N+3 . 


-2acr 


From Lemma 6.4 and (15.481) . we have 

^I|l 2 (d) 


llViL " 2 


IIVxfltH? 


< c. 


i 2 (£>) 


(5.48) 


(5.49) 


Combining Lemma 6.4, (15.4711 and (15.491) we obtain 

iivx1ni 2 !L - 


II L 2 (D) 


II L 2 (D) 


|V X Hi 112 


L 2 (D) 


+CrS r 


II ViL , 112 


L 2 (D) 


IIV X Fftllr 


L 2 (D) 


< C{(1 -r)<5 _(1_r) 1 + rS r 'jO^" 2 ) + Cr5 r \ 
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We now choose (|,2), combining (15.461) . (15.471) and (15.491) we have 


l|Vx^||l 2(D) 


^ \mh {D) \m\ 2 Hr (D) ii v x H t \\ 2 LP ( D) 

Cl ||Vxi7 t ||| 2(D) C2 \\VxHt\\l 2{D) C3 ||Vx^||| 2(d) 

> C — c\{{\ — r)5~ < ' 1 ~ r ' > +rS r }0(t~ 2 ) — CrS r — c^t 1 ~p 

> C — cirb r , - < r < 1, r » 1. 


Hence from Lemma 6.6, we have 


liminf \I p (t, h B {p))\ > c > 0. 


6 Appendix 

6.1 Construction of the oscillating-decaying solutions A 
and B 

In this subsection, we show how the scheme in m can be used to derive the 
oscillating-decaying solutions A and B. Recall that E and H satisfy equation 
& therefore we need to derive estimates of the higher derivatives for A and B. 

Note that the main term of w A t b t N u (resp. w^ t b t N A is xt{x , )Qte lTX '^e~ T ^ x ' u ~ t ' ,A * f ' x '' ) b 

(resp. Xt(A)Qte lTX ^er T<x ' UJ ^ t ^ A A x ''>lj), which can be directly differentiated term 
by term since it is a multiplication of smooth functions. So we can calculate E 
and H directly. For convenience, we denote w = w Xttb! t,N,u 7 = 7**, b,t,N,u{x, t). 

Without loss of generality, we can use the change of coordinates to assume t = 0, 
w = (0,0,1) and r] = (1, 0,0), £ = (0,1, 0). Define 

Qa := e"*™'*' L A {e iTX '<'•), Qb := e~ iTX '<' L B {^ rx '<' ■) 

where x' = (xi,X 2 ), = (£1,^2) with |£'| = 1 and La,L b have been defined 

by (13.51) and (13.61) . In the following, we will give all the details for the higher 
derivatives of E and H. 

In [15] , the authors used the phase plane method to get a first order ODE 
system and we want to decouple the equation in order to solve it by direct 
calculations. The method of construction the oscillating-decaying solution is 
decomposed into several steps: 

Step 1. As mentioned before, we set <5,4 = e~ lTX ''^LA(e lTX ' '£'■), Q b := 
e~ lTX L B {e lTX •) and solve Qava = 0, Q b v b = 0. In the following calcu¬ 
lations, we only need to consideiMjA'iM = 0 since Q B v B = 0 will follow the 
same calculations. Let Qa = C a Qa be the operator which satisfies the leading 
coefficient of d | is 1 and the existence of Ca is given by the strong ellipticity of 
La and we need to solve Qava = 0 (the same reason for the operator Q B and 
Q b ). Now, We introduce the concept of the order in the following manner. We 
consider r, 83 are of order 1, di, 82 are of order 0 and X 3 is of order —1. 

Step 2. Use the Taylor expansion with respect to X 3 , we have 

x N ~ 1 

Qa(x',X3) = Q a {x', 0) + • • • + ^ 3 _ -Q! 83 1 Qa(x' } 0) + R 

= Q\ + Q a + • • • + Q a N+1 + ft 
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where ord(Q^) = j and ord(i?) 


— N. Since we hope that Qava = 0, we have 


Q\va — ~(Qa + Qa + • • • + Qa N+1 + R) v A •— f■ 


Step 3. Following the paper na, we denote D 3 = —id 3 , p = (£ 1 ,^ 2 ,0) and 
(a, b) = ({a,b) ik ) for a = (ai,a 2 ,a 3 ) and b = (bi,b 2 ,b 3 ), where ( a,b) ik = 
C^ kl ajbi with C-j kl being the leading coefficient of the second order strongly 


elliptic operator La- If we set W 


Wl 

w 2 


where 


I Wi = v A 

\^ 2 = ~ T ~ l ( e 3 ; e 3 ) X3=0 D 3 v a - (e 3 , p) x 3=0 v A 
and use / = — (Q\ + Q\ + • • • + Q^ N+1 + R)va, then W will satisfy 


D 3 W 


tK a W + 


0 

r- 1 (e 3 ,e 3 ) X8=0 / 


(■ tK a +K a + --- + K a n + S)W 


where K A is a matrix in depending of x 3 which can be diagonlizable by the 
property of the strong ellipticity of La- Note that each K A, s only involves the 
x' derivatives with ord(A'' 4 ) = j, ord(S') = — N — 1. It is worth to mention that 
with the help of such special W, then we can solve the ODE system explicitly. 
Step 4. Decompose K A such that 


K A = Q- x K a Q 


K A JL 
0 K a 


where spec(K A ) C C± := {±Im\ > 0} (the existence of K A and Q were showed 
in [T5])- If we set W = Q _ 1 W, then 


D 3 W = ( tK a + K 0 + --- + K_ n + S)W, 


Step 5. If we write W = (I + x 3 A 1°) + with A®, being differ¬ 

ential operators in d x > (their coefficients independent of x 3 ), then 

D 3 W {0) = {tK a + {K 0 -txsAWra + txsK^AW - B^K* 

+K A B {0) + L4 (0) ) + K>_ 1 + • • • }W (0) 

:= ( tK A + K 0 + K'_ 1 + • • ■ )W (0) 

where ord(AT'_i) = —1 and the remainders are at most —2. We choose A( 0 \B( 0 ' > 
to be suitable operators and use the same calculations in m, then we will get 


K 0 = 


tfo(l.l) _ 0 
0 K 0 (2,2) 


36 










to be a diagonal form (here we omit all the details). 

Step 6. Finally, following step 5, we can write 

W = {I+ x 3 A {0) +t~ 1 B < ' 0) )(I+ x 2 3 A {1) +t- 1 x 3 B {1) 

x (J + x% +1 AW + + T-txg-'cWjwW 

with suitable A^\B^ and for j = 0,1,2, ■ ■ • , N (C*- 0 -* = 0), then W ^ 

satisfies _ _ 

D 3 W (n) = {tK a + K 0 + • • • + K_ n + 5}VF (Ar) , 

with all K_j are decoupled for 0 < j < N and ord(S') = —N — 1 . If we omit 
the term S , we can find an approximated solution of the form 

JV+l 

< 4 "’ = E 

3=0 


satisfying 


Atfi * 0 = + K 0 (l, !) + ••• + K_ jv(1, 1)K 


(JV) 


and each 0 ^^ has to satisfy 


d 3 v { o N a = tK&W > 

d 3 ^-i] A = tK +v < L^ a + i^o(l, 1 )Dq|a, 


ft 


V 


W| 

0,A 1*3 

W , 

-i,aI*3=o 


=0 = xt{x')b, 
= 0, 


- tK aAN) ,y,N £ o n-w -w , 0 

— ril + t ; _jV-l,A w Zjj=0 1 1"-J,A> v -JV-1,aI* 3=0 — u ; 

where Xt(x') £ C(j > °(]R 2 ) and b £ C 3 . Thus, by solving this ODE system we can 
get the following estimates: 


£3 5 “'(^,a)IIl 2 (R3.) < CT 


-P-j-1/2 


( 6 . 1 ) 


for 0 < j < N + 1. Moreover, if we set = VA 
-J - ’ A 0 


-AN) 


then it satisfies 


'vf° - {tK a + K 0 + --- + K_ N }V[ N) = R, 


T AN) 

V A 1*3=0 


Xt{x')b 

0 


where 


II^IIl2(r 3 ) < cr N 3/2 . 


Step 7. Finally, if we define the function v A = 


Vi 

Vi 

v 3 


, with Vj being the jth 


component of the vector Q(I + x 3 A^ + r 1 B^)(I + x^A^ + r 1 x 3 B ( - 1 ' ) + 
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T- 2 C«) •••(/ + X» +1 AW + T -1 x N B (N) + T -2 x N-l c ( N)) y{N) &nd get ^ = 

exp(zr*' • £')va, we will get that 


wa = Q exp(zrx / ■ £') exy{iTX 3 K^:{x'))xt{x')b + exp(irx / • £ 7 )r (x, r) 

= Qexp(irx' ■ £') exp(-iTX 3 (-Kf(x')))xt(x')b + r(x, r) 

and 

wa\x 3 =o = exp^ra/ • £'){xt(x')Qb + Po(x r , r), 

where /3o(x',t) = r(x',0, t) is supported in supp(xt). Note that the function 7 
comes from the combination of u^^’s, for j = 1.2,--- , TV + 1. Now, we derive 
higher derivative estimates for the oscillating-decaying solutions, back to see all 
the i)^ A ’a separately. In fact, only need to see From the estimate (16.11) . 

we know that the estimate is independent of the derivative of x' variables, all 
we need to concern is the 83 derivative. From the equation 

Dsv^a = rKfv^A + K 0 (l, 1 )*>£2 ( 6 . 2 ) 

and the standard regularity theory of ODEs(ordinary differential equations), we 
know that £ C°° if all the coefficients are smooth. Moreover, note that 

K + independent of x 3 , then we can differentiate (IQ) directly, to get 

DIv^a = D 3 [t^vW a +K 0 (1,1 )vI”2 ] 

= + (D 3 K 0 (1 ,1 ))v ( 0 N 2 + K 0 (l, 1 W™ 

= T^fvW A + rKfK 0 (l, 1 )f)g + (D 3 K 0 (1, 1 ))fig 

+TK 0 (l,l)Kfv { 0 N 2. 

Thus, we can obtain that 

ll a: 3 a “' 9 3( fi - J ]L 1 ) A)IU 2 (K^) < cT- 0+ri ~ 3/2 , 


for all rf < 2. Inductively, we have 


ll^3 a “' a 3(^- 7 I?A)llL 2 (K3) < 


-P+1 7-3/2 


for all 77 £ N. Similarly, for other v^\ w ifh 2 < j < N + 1, we can get similar 
estimate in the following: 


x 3 d x' d 3(^A)\\^m) < ct 


n-p-j- 1/2 


Vr] £ N U {0}. Therefore, T satisfies 

II^FlU^n.) < cr H-3/2 e -T (s -i)A 


on fl s := {x 3 > s} D fi for s > 0 and V|a| £ N U {0}. Note that since each 
v_p A ’s are smooth, we can get the smoothness of R and 

I|S“ j R|U 2 (R 3_ ) < crM-^- 3 / 2 
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for all |cc| eNU{0}. Furthermore, we have that 

\\dZ{QAVA)\\L*(n 0 ) < ct^~ n ~ 1/2 . 

Step 8 . Now let u = w + r = e lTX v + r and r be the solution to the boundary 
value problem 

L A r = -e lTX '^' Qava in H 0 

r = 0 on OCIq 

However, note that Ho = { 1 E 3 > 0} l~l Cl is not a smooth domain since dClo = 
({X 3 = 0} n Cl) U ({X 3 > 0} l~l dCl). Note that the oscillating-decaying solution 
exists in the half space, from the construction, we know that the solution is 
independent of the domain Cl. Let Cl C R(j_ be a open bounded smooth domain 
containing Cl with {X 3 = 0} n Cl C dCl, from the construction, it is easy to see 
the form of oscillating-decaying solution does not depend on the domain Cl, then 
we can extend r to be defined on Cl and call it f(x). Here we can also extend 
va to be defined on H, still denote va and all the decaying estimates will hold 
since our estimates were considered in R+, then we have 

L A f = —e tTX '^' Q a va in Cl, 

r = 0 on dCl. 

Note that all the coefficients are smooth, we apply a well-known elliptic regu¬ 
larity theorem (Theorem2.3, |T|), then we will get f £ C k (Cl) Vfc (recall that 
dCl £ C°°) and 

ll^lliP' + 1 (f1;R 3 ) < c||<3j4i;A||fJ*:(n;]R3). 

Hence \\d^r\\ L 2 ^ 0 ^ < ||<9“f < crl“l _Ar+1 / 2 for all |a| < k, Vfc £ N. Simi¬ 

larly, we can construct the oscillating decaying solution for LbB = 0. Then we 
represent A and B to be two oscillating-decaying solution in the following form: 

W Xt,b,t,N,u> + r xt,b,t,N,u: in 

e iTX <{xt(x')Qt(x')b + t t bNuj } on E t (w), 

W Xt,b,t,N, U + r xt,b,t,N,u in n t(u), 

e iTX< {xt{x')Q t {x')b + t J on E t (w), 


Xt(x , )Q t e iTX <e- x ^- t ^ A ^b + l^ N Jx, r), 
Xt(.x')Q t e ixx <e-< x “-* A ?t x ')b + 7^,&,*,*,„(*, r), 

7 £ u b,t,N,u and ixt,b,t,N,w satisf y (EH) and EH- 

6.2 Well-posedness and L p estimate for the anisotropic 
Maxwell system 

In the following, we would list the eigenvalue property and well-posedness results 
of the following problem: let Cl C R 3 and K <g Cl, 


where 


W Xt,b,t,N,u 
B 


w 


Xt ,b,t,N,u 
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(6.3) 


V x E = ik/iH in f l \ K 

V x H = -ikeE + J in Q \ K 

v x E = / on dfl 

v x H = g on dK , 


where p,, e are symmetric and positive definite matrix-valued functions. More 
precisely, we assume there exist constants /io, /ri, Ao, Ao > 0 such that 

[A 0 I < e(x) < A 0 I. 

These well-posedness for the isotropic Maxwell systems can be found in Theorem 
4.18 and 4.19 of jlB]. However, we have the same result under our assumption 
(16.411 following the arguments in [13]. Let 

X = £ H(curl;fl\ K)\v x u = 0 on dfl and ut £ L 2 ( dK ) 3 on <9 At| . 

Definition 6.1. We say (E, H) or If is a weak solution of (16.31) if E £ X and 
satisfies 

V x E, V x &) n \ K -k 2 ( 7 E, (f>) U \ K = {ikJ , <fa) aK - W> £ X, 

(6.5) 

and v x E = f on dfl, where (f>r = {v x (j>) x v and (•, •) denotes the standard 
Hermitian inner product of L 2 space. Moreover, if (16.51) fails to have a unique 
solution, then k is called an eigenvalue or a resonance of m- 

Lemma 6.2. There is an infinite discrete set E of eigenvalue kj > 0, j = 
1,2,... and corresponding eigenfunctions Ej £ Ho(curl;Q), Ej 0, such that 
116.50 holds with J = 0 and f = g = 0 is satisfied. 

From the above lemma, we have the following theorem. 

Theorem 6.3. For k E, there exists a unique weak solution ( E , H) £ 
H(curl;n\K) x H{curl]kl\K) of \6.3\ ) given any f £ H~ 1 ^ 2 (Diw 1 dfl), g £ 
H~ 1 ^ 2 (Diw,dK) and J £ H~ 1 (fl\K). The solution satisfies 

\\m L 2 (n\K) + W H IU 2 (n\if) < C{\\f\\ H -i/ 2 ( Div . dn )+\\g\\ H ^/ 2 ( Div . dK) + \\J\\ H - 1 ^^) 
for some constant C > 0, where 

H- 1 ' 2 C Div, T) := {/ £ H- 1 / 2 (r) 3 | w ■ / = 0, Van • / £ H~ 112 (T) } , 

T = dfl or dK. 

In the following, we state the L p theory for the anisotropic Maxwell’s system. 
For this purpose, we define a bilinear form 


B a {E , F) := / (A(x)V x E{x)) ■ (V x F(x))dx + M E{x) ■ F{x)dx 


in 
1 1 


for all E £ HQ' q (curl,fl) and F £ ( curl,fl ) with —I—- = 1. We only 

state L p estimate in the following theorem, but we do not prove the theorem. 
For more details, we refer readers to read [7]. 
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Theorem 6.4. Jffl] Let LI be a a smooth domain. Suppose that A = A(x) is 
a real symmetric matrix with smooth entries and satisfies the uniform elliptic 
condition 

A |(| 2 < ;4(a:)f • f < A|(| 2 , for all £ G R 3 , 

for some constants 0 < A < A < oo. Assume q is some number satisfying 
2 < q < oo. Under the condition 


inf sup \BaE,F)\ > — > 0 
ll F lll,?' =1 ll^lll.O = l K 


the Maxwell’s systems of the equations 


V x (AV xE) + E = Vxf + g 

is uniquely solvable in iJg ’ 9 (curl, LI) for each g G L q (LI) and f G L q (LI) and 
the weak solution satisfies 


\\E\\L*fn) + IIV x E\\ Lql{n) < K{\\f\\ Lq , m + || 5 ||^ (n) }, 

where K is a positive constant depending on p. 

We end up this appendix with the following lemma on the embedding related 
to the Sobolev-Besov spaces, for more details, see 1111- 

Lemma 6.5. Let u G L P (D) such that V • u G L P (D) and V x u G L P (D). 
If v x u G L p (dD), then also v ■ u G L p (dD) for p G (l,oo). If in addition 
1 < p < 2, then u G B P £ 2 (D) and we have the estimate 

V 

IMIb^CD) — C{IMIi p (I}) + ll CMr HI LP(D) + ||V • m||lp(_d) + \W x M |l LP(dD)} 

P 

where the Sobolev-Besov space B p ’ q (D) := [L p (D),W l ’ p (D)\ a ^ is obtained by 
real interpolation for 1 < p, q < oo and 0 < a < 1 . 
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